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ABSTRACT 

We use the methods of group theory to reduce the equations of motion of two 
spin systems in (2+1) dimensions to sets of coupled ordinary differential equa- 
tions. We present solutions of some classes of these sets and discuss their physical 
significance. 



Lcs methodes de la theorie des groupes sont utilisees pour reduire les equations 
du mouvement de deux systcmes de spins de dimensions (2+1) a des systemes 
d'equations differentielles ordinaires. Les solutions de certaines classes de ces 
systemes sont presentees et les aspects physics sont discutes. 



1. Introduction 

In this paper we look for solutions of the equations of the Landau-Lifshitz 
model (with, perhaps, nonvanishing anisotropy) and of a nonlinear vector diffusion 
equation. The equations are given, respectively, by 



f^.-x# (1.1) 

and 

— * 

^ = F-0(0.F), (1.2) 

where F is given by 

F^A$+{Acj>3 + B)e3, (1.3) 

where 63 is a unit vector in the 3'^'^ direction in the space and satisfies 0.0 = 1. 
A and B are possible anisotropy coefficients. 

The motivation for this work comes from the original observation made by 
Landau and Lifshitz'^' in their study of the ferromagnetic continuum. They pointed 
out that for phenomena for which substantial spatial variations occur only over a 
large number of lattice spacings, we can use the continuum approximation. They 

showed that a ferromagnetic medium is characterised by the magnetisation vector 
M (like the vector <p above) which precesses around the effective magnetic field 
and so obeys, what is now called the Landau-Lifshitz equation, namely (1.1). Since 
the original work of Landau and Lifshitz many papers have been written on the 

r2] 

subject and the equation has been modified by the inclusion of various additional 
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terms to F. It has been used to describe the dynamics of magnetic bubbles in a 
ferromagnetic continuum and also of vortices in Hell or in a superconductor . 
Various studies of the dynamics of such topological soliton-like structures have 
been performed both theoretically and experimentally'^' and they have exhibited 
many interesting, and perhaps unexpected, phenomena - like the skew deflection 
of these structures under the influence of a magnetic fleld gradient which resembles 
the more familar Hall motion of electrons in external magnetic and electric flelds . 

A recent work of Papanicolaou and Tomaras , as well as some earlier work of 
other people'^' has shown that many experimentally observed facts can indeed be 
explained using the Landau- Lifshitz equation. Much of the work involved deriving 
various conserved quantities describing these structures and then using them to 
restrict the description of the dynamics. All this work has provided further evidence 
as to the relevance of the Landau-Lifshitz equation to the description of physical 
phenomena. However, as the Landau-Lifshitz equation is quite complicated, only 
some results were obtained in an analytical form. Most more recent studies 
involved numerical simulations. 

The vector nonlinear diffusion equation (1.2) has less obvious physical appli- 

[^] 

cations but it has been used in the study of phase ordering kinetics where one 
investigates the time evolution of a system quenched from the disordered into an 

rgj 

ordered phase. This topic has attracted considerable attention in recent years . In 
fact, it has been shown that many features of phase ordering in systems supporting 
topologically stable defects (for example, in systems described by the O(A^) vector 
model in d dimensions with d < A^'^°' or in two and three-dimensional nematic 

[12] 

liquid crystals can be understood theoretically by investigating the dynamics 
of the numerous topological defects generated during the quench. A special and 
interesting case is that of the 0(3) model system in 2 spatial dimensions. It sup- 
ports topologically stable, but non-singular objects which, in the condensed matter 
community language, are called topological textures. Such systems were studied 
numerically in ref. [8]. 
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Given the paucity of analytical results for both equations (1.1) and (1.2) (es- 
pecially involving the dynamics) one of the aims of this paper is to see what time 
dependent solutions can be found using the group theoretical method of symmetry 
reduction . This method exploits the symmetry of the original equations 
to find solutions invariant under some subgroup (the classic example one can give 
here involves seeking solutions in three dimensions which are rotationally invari- 
ant). The method puts all such attempts on a unified footing and it has been 
applied with success to many equations . The method gives equations whose 
solutions represent specific solutions of the full equations; the solutions are deter- 
mined locally and the method does not tell us whether these solutions are stable 
or not with respect to any perturbations. 

In a recent paper '^'^'j two of us (PW and WJZ) together with M. Grundland, 
have apphed this technique to looking for solutions of the relativistic CP^ model. 

In this paper we investigate solutions of (1.1) and (1.2). We are particularly 
interested in time dependent solutions; all time independent solutions of (1.1) and 
(1.2) (when there is no anisotropy) are also the time independent solutions of the 
relativistc model and so can be found in ref [17]. 

Like in the relativistic CP^ model studied before, in order to perform the 
symmetry reductions, we have to decide what variables to use. To avoid having to 

use the constrained variables (0) it is convenient to use the W formulation of the 
model which involves the stereographic projection of the sphere 0-0=1 onto the 
complex plane. In this formulation instead of using the fields, we express all the 
dependence on in terms of their stereographic projection onto the complex plane 

— * 

W. The fields are then related to W by 

w + w* , w-w* , i-\wf , 

01 = \ — ;2' ^2 = I ■ — ;2' ^3 = \ — -n- (1.4) 

i + iiyr i + ii^r i + |w^r 

To perform our analysis it is convenient to use the polar version of the W 
variables; i.e. to put W — R exp iQ and then study the equations for R and Q. The 
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advantage of this approach is that the equations become simple; the disadvantage 
comes from having to pay attention that R is real and Q should be periodic with a 
period of 27r. (If the period is not 27r then the solution may become multi-valued) 
Thus if we find solutions that do not obey these restrictions, then these solutions, 
however interesting they may be, cannot be treated as solutions of the original 
model. 

In the case of the Landau- Lifshitz equation the equations for R and Q take the 
form 

dtR - 2 |j~^^| {dxQdxR + dyQdyR^ - R{dxxQ + dyyQ) = (1.5) 



and 



dtQ ^B + A 



1 — R^ dxxR + dyyR 



l + i?2 R 



(1.6) 



while for the diffusion case they are respectively 

dtQ - 2 ^[^~^^2) (dxQdxR + dyQdyR) - {dxxQ + dyyQ) = (1.7) 



and 



dtR +BR + AR——^-dxxR - dyyR 



1 + R^' 

[idM + (dyQf) {{dxRf + (dyRf) = 0. 



(1.8) 



Note, that, in the Landau-Lifshitz case, if we put R = 1 the equations become 
AQ = and dtQ = B which have a very simple solution, and in the diffusion case, 
we have to set 5 = and then we end up with dtQ — AQ = as the equation for Q. 
The latter case is the nonrelativistic analogue of what was found in the relativistic 
case where R — 1 reduced the equation for Q to the linear wave equation for the 
phase Q. 
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In the next section we determine the symmetry group of our equations (1.6) 
, (1.5) and of (1.7) and (1.8) . In the following sections we solve the derived 
equations and discuss their solutions. 



2. The Symmetry Group and its Two Dimensional Subgroups 

The symmetry group of our systems of equations, respectively (1.6) and (1.5) 
and (1.7) and (1.8) , can be calculated using the standard methods '""""^^"^'^^ . 
We actually made use of a MACSYMA package that provides a simpUfied and 
partially solved set of determining equations. 

Solving the determining equations we find that three different cases must be 
distinguished: 

1. A — B — 0, i.e. the anisotropy is absent. The Landau-Lifshitz equation 
and the diffusion equation have isomorphic symmetry groups, consisting of 
translations in space and time directions, rotations in the x, y plane, dilations 
and a group of 0(3) rotations between the components of the field 0. The 
corresponding Lie algebra Li has the structure of a direct sum 



Bases for these two algebras are given by the following vector fields, acting on 
space-time and on the fields in the {R, Q} reahsation of eq. (1.5)-(1.8): 



Li = s(2,l) e 0(3). 



(2.1) 



s(2,l) 



Po = dt, Pi = d, 



Xl 



P2 = d, 



(2.2) 



L = -xdy + ydx, D = 2tdt + xdx + yd^ 



su{2), 



Y 



X 



Z 




(2.3) 



2. A^Q. 
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The symmetry algebra for both equations is reduced to 



L2^{Po,Pi,P2,L} ® {Z}, 



(2.4) 



i.e. the dilations are absent and the only rotations left are those around the 
third axis {i.e. around (ps). 

3. A^O, S ^ 0. 

The symmetry algebra for the dissipative equations (1.8) and (1.7) is still L2, 
as in (2.4). That of the Landau- Lifshitz equation is 



In order to perform symmetry reduction we need to classify the subalgebras 
of the symmetry algebras Li, L2 and L3. We wish to reduce equations (1.5)- (1.8) 
to ordinary differential equations. To do this, we will require that the solutions 
are invariant under a two-dimensional subgroup of the symmetry group. In order 
to do this systematically we need to derive a classification of the two dimensional 
subalgebras of the symmetry algebra. Moreover, we can restrict ourselves to sub- 
algebras, all elements of which act nontrivially on space-time, i.e. which do not 
contain any rotations in 4> space. 

The subalgebra classification can be done in an algorithmic way'"'; the results 
are quite simple and we present them without a proof. 

1. A — B — 0. Every two-dimensional subalgebra of Li, each element of which 
acts nontrivially on space-time, is conjugate under the action of the group of 



L3^{Po,Pi,P2,L,D} ® {X,Y,Z}, 



(2.5) 



with 



D = 2tdt + xdx + ydy + 2BtdQ. 
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inner automorphisms to one of the following ones 



^2,1 
^2,2 

^2,4 

^2,5 
^2,6 

M,i 

^2,8 



{Pi + aZ,P2 + hZ}, 
{L + aZ,Po + bZ}, 
{Po + aZ, P2 + bZ}, 



{Pq - vPi + aZ, Pa + bZ}, Vy^O, 

{D + bL + aZ,Po}, by^O, 
{D + aZ,PQ}, 



(2.7) 



{D + aZ,L + bZ}, 
{D + aZ, Pa}- 



The parameters a, b and v are arbitrary real numbers. In some cases their 
ranges can be further constrained but that is not important for our purposes. 



Every two-dimensional subalgebra of the considered type is conjugate to one 
listed above as Aa,!, ...Aa,4- 

3. A^O, S ^ 0. 

For the dissipativc equations (1.8) and (1.7) the subalgebra classes are repre- 
sented by Aa^i, ...A2A- Foi' the Landau-Lifshitz equations (1.5) and (1.6) they are 
represented by ^a,i, •••^2,8 with D replaced by D and Pq replaced by Pq = Po + bZ. 

We can now proceed to perform various reductions. We are particularly inter- 
ested in reductions that do not result in time independence as these were already 
studied in ref [17]. 



2. Ay^O. 
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3. Solutions of the Landau-Lifshitz Equation 



3.1. General Procedure 

Our aim is to solve the Landau-Lifshitz equations (1.5) and (1.6) , using the 
method of symmetry reduction. This involves assuming that a solution is invariant 
under a subgroup Go of the symmetry group G, namely one of the two dimensional 
groups corresponding to the algebras ^2,1, • • • , ^2,8 of (2.7). The assumption makes 
it possible to reduce the partial differential equations (1.5) and (1.6) to a pair of 
coupled ordinary differential equations. Whenever possible, we decouple them and 
find exphcit solutions for the functions R and Q, hence for W, and finally for the 

— * 

vector figuring in (1.1) . 

For all 8 algebras in (2.7) the invariant solution will have the form 



where ^ and (5 are explicitly given and R{^) and a{S,) satisfied coupled ordinary 
differential equations obtained by substituing (3.1) into (1.5) and (1.6) . 

The reduced equation (1.5) is 



(V0^^^C+[2 ^ J (VO^i?^ + Ae]a^ = -|6-2-^-^i?5(VC,V/?)-A/?. 



R{x,y,t)^R{0, 



Q{x,y,t) 



a{0 + P{x,y,t) 



(3.1) 



(3.2) 



For algebra ^42,1 we have 



we = AC = A/3 = 0, 



(3.3) 



and so (3.2) reduces to R^ — 0. 

In all other cases we have {Vif + 0. Eq. (3.2) is a first order linear inhomo- 
geneous equation for a^. We can integrate it explicitly and obtain in terms of 
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R, whenever the functions ^ and /3 satisfy 



d f h 



where 



m = 1 



h 



(ve,V/g) ' 



hA(3 



(3.4) 



for = 0, 
for 7^ 0. 



(3.5) 



Conditions (3.4) are always satisfied for the algebras ^2,2 ■ ■ ■ ^2,6, not however 
for ^2,7 and ^2,8- When conditions (3.4) are satisfied, we can integrate eq. (3.2) 
once to obtain 



sd + ii 



2\2 



1 + 



h R'' ^ i?2 

where S is an arbitrary real integration constant and where we have 



(3.6) 



V — 

V - 



62 + 1' 



— 1+^2, 
i/ = 0, 



for 
for 
for 



^2,4 
^2,5 
^2,6 



(3.7) 



for A2,2,>l2,3- 



Equation (1.6) for algebras ^2,2, • • • ,^2.8 is reduced to a second order differ- 
ential equation for i?(4), that also involves a^(C)- For reductions corresponding to 
Lie algebras ^2,2, ■ ■ • , ^2,6 we can substitute from (3.6) , to obtain an ordinary 
differential equation for i?(^) alone. To transform this equation to a standard form 
we put 



^(0 



(3.^ 



The equation for U (77) is then written as 
1 



, 1 



Equation (3.9) can be integrated in terms of elliptic functions if p, q and m are 
constants. This is always the case for algebras ^2,3, ■ ■ ■ , ^2,6- In the case of algebra 
^2,2 this is true if we set ^4 = 0, S = h. 



^^'(C/+l)(C/-l)3+p^^^^+5C/+m(C/-l)2. (3.9) 
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Eq. (3.9) has a first integral that we can write as 

= -4S^U^ + KiU^ + KU^ + K2U + (3.10) 

where K is an integration constant, and the constants Ki,K2 and Ks are related 
to the coefficients S,p, q and m in (3.9) . 

In this article we restrict ourselves to solutions of the Landau-Lifshitz equation 
that are obtained by solving (3.10) . 

We shall first discuss solutions of (3.10) in general, then run through algebras 
^2,2, • • • , ^2,6 and specify the values of the coefficients in (3.10) in each 
well as the independent variable rj. 

Algebra ^2,1 leading to a first order equation, will be treated separately. 
3.2. Solutions of the elliptic function equation 

We shall call (3.10) the "eUiptic function equation". Its solutions are of course 

[19] 

well known . We shall however list those that are relevant in the context of 
solving (3.9) , and more importantly, the Landau-Lifshitz equation. 

Several comments are in order here: 

1. The functions R{ri) must be real (and nonnegative) , hence U{rj) must be real 
and nonnpositive. 

2. For 5* 7^ the coefficient of the highest power of U in (3.10) is nonnegative. 
this means that all real solutions of (3.10) are nonsingular. 

3. For S — 0, Ki in (3.10) the real solutions of (3.10) can be singular. Since 
we are really interested in the fields 0i we note that singular solutions of U 
will give regular functions 

4. In general, equation (3.10) is solved in terms of Jacobi elliptic functions. 
However, these reduce to elementary functions whenever the polynomial on 
the right hand side has multiple roots, or when S — Ki — 0. 



11 



Let us run through individual cases. 

I. Sy^O 

We rewrite (3.10) as 



C/2 = -iS^U - Ui){U - U2){U - Us){U - U4) (3.11) 

1. Ui<U <U2 = U3 = U4<0 

this is an algebraic solitary wave, equal to U2 for 77 — > ±00, and dipping 
down to Ui for 77 = 770. 

2. Ui = U2 = Uz<U <Ua<Q 

Also an algebraic solitary wave, rising io U — U^^ for 77 = 770, equal to Ui 
for 77 — > ±00. 

3. Ui<U <U2^Ui< U4, U2<0 

Ui,)^U2 iU,-U2)iU2-U,) 



{Ui Ih) cosh^ fi{r] - 770) - (C/2 - Ui) (3.14) 



/i=>Sv/(C/4-[/2)(t/2-[/l) 

4. C/i < C/2 = C/3 < [/ < C/4 < 

C/3 + '^^-^'"^^-^^> (3.15) 

with /X as in (3.14) . 
The last two solutions are solitons, the first one a well, the second a bump. 
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5. Ui<U <U2<Uz = U4, U2<Q 

(C/4-C/2)(C/4-C/l) 



U{rj) =U4 



{U2 - Ui) sin^ n{r] - rjo) + U4-U2 (3.16) 



^x^S^/{U4-U2){U4-Ul) 
6. C/i = C/2 < C/3 < [/ < C/4 < 



(C/4 - Us) sin' i^{ri - rio) + U3-U1 (3.17) 
pL^S^/(U4-Ul)(Us-Ul) 



7. Ui<U <U2<U3< C/4, U2<0 

{U4-U2){U4-Ui) 



U(v) =C^4 



{U2 - Ui)sn^{ii{r] - 770), k) + U4- U2 



(3.18) 



8. Ui<U2<U3<U <U4<0 

""^^^ = + (c/4-c/3).nM^-^o) + c/3-c/i ^'-''^ 

with A;^ and as in (3.18) . 

9. Ui <U <U2 <0, C/3,4 = p±tq,q>0 

, . _ {MUi - NU2)cn(ii(v - 770), k) + MUi + TV^/^ 

(M - iV)m(/i(r7 - ?7o), k) + Af + 
M^^{U2-pf + q\ N^^{Ui-pf + q^ (3.20) 

^^^^ 

AMN 

Solutions (3.16) , . . . , (3.20) are periodic. All the elementary solutions can be 
viewed as limits of solutions (3.18) , (3.19) and (3.20) . 
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//. S = 0, Kly^O 

Set 



1 



/x=-V|i^i|([/3-[/i (3.21) 



1. < 0, C/i = C/2 < [/ < C/3 < 

[/ = C/g - (C/g - U2) tanh^ ij{r) - 770) (3.22) 

2. i^i < 0, t/ < [/i = [/2 < t/3, t/i < 



(^3 - Ul) 

tanh^ /x(?7 — tjq) 



U-Us- ' . (3-23) 



3. Ki<Q,U <Ui = U2 = U2,<Q 



4. < 0, [/ < C/i < C/2 = C/3, C/i < 



^3-^1 
sin^/x(?7 - ?7o) 



^ = ^3 - . ;\ . (3.25) 



5. < 0, C/i < C/2 < [/ < C/3 < 



U^Us-iUs- U2)sn\i^{ri - 770), A;), A;^ = (3.26) 

U3 - Ul 



6. /sTi < 0, C/ < C/i < C/2 < Us, Ui<0 



U = Us- J^^ ^\ (3.27) 
sn'^^li^r] - rjo),k) 



k as in (3.26) 
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7. Ki > 0,Ui < U < U2 = Us = 

1 



cosh iJi{ri — Tjo) 

8. i^i > 0, t/i < [/ < C/2 < < t/3 

U^{U2- Ui)sn^{f,{ri - 770), k), = 

9. Ki <0,Ui< 0, C/2,3 = p ± ig, g > 

U^Ui+A 



1 - m(/x(77 - r]o),k) 



///. 5 = 0, /Ti = 0, 7^ 
1. X > 0, [/ < c/i < < C/2 



2 



U^Ui-(U2- Ui) sinh" - ^0) 



2. X > 0, [/ < C/i = C/2 = 

[/ = - exp(-V^(r7 - ?7o)) 



3. X < 0, C/i < [/ < C/2 < 

u^Ui + {U2- Ui) cos' ^^{v - m) 

IV. S ^ Ki^ K ^0,K2^0 

U = -^^ + ^iv-Vof, K2<0, Ks<0 
V. S = Ki = K = K2 = 
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3.3. Individual reductions 



1. Algebra /l2,i. 

This is an exceptional case when (3.2) implies — 0. We find that the only 
solution for W of (1.4) is 

W^Rqc'^, Q^ax + by+(B + ^—^{a^ + b^ + A))t + aQ (3.36) 

1 + Rq 

where Rq and ao are integration constants. 

2. Algebra ^2,2- 
We find 

W ^ R{p)expi[a{p) + a(t) + bt], ^ = p (3.37) 

where p and are polar coordinates. The singlevaluedness of W requires a to be 
an integer. The phase a{p) and variable rj satisfy 

(1 + -R^)^ 

(^p{p) = S — — ' ^ = lnp (3.38) 

(see (3.6) ). For the function U{r]) of (3.8) we obtain the elliptic function equation 
if and only if we set 

A = 0, b^B (3.39) 

{A and B are defined in (1.6) ). 
We have 

in (3.10) . 
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For S* 7^ cq. (3.40) implies that wc can have two negative and two positive 
roots in eq. (3.11) or two negative roots and two complex conjugate ones. These 
cases lead to real solutions, namely (3.16) , (3.18) and (3.20) . Note that all of 
them are periodic. In particular, for a = eq. (3.11) always has a double root 
Us — U4 — 1 and reduces to 

C/2 = -AS'[U^ + (i + JL)u + l] (3.41) 



For S — 0, K Aa^ we obtain the equation 



C/2 = 2(a2 - ^)U{U - Ui)(U - U2), 
U1U2 = 1, t/i + t/2 = 



X - 4a2 ■ 



(3.42) 



The relevant solutions of (3.42) in this case are: 

1. 2a'^ <K < 4a2, Ui < U2 < 

solution (3.29) (with U3 = 0). For K = 20^ we have U = -1. 

2. K> 4a2,f/i <U2<0 

Solutions (3.26) , (3.27) , (3.22) and (3.23) (all with t/s = 0). 

3. K> 4a2, < C/i < C/2 

Solutions (3.26) (with U3 U2, U2 Ui, C/i = ). 

4. K > 4a2, Ui^2 = p±iq,q> solution (3.30) . 

For S — 0,K — Aa^ (3.10) reduces to an elementary one and its solution is 



where Rq is an integration constant. 
3. Algebra ^2,3 and ^2,4 
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The reduction formulas in both of these cases are 



W ^R{Oexpi[a{0-at-by], C^x + vt, rj ^ C (3-44) 

with V — and ^; 7^ for the algebras ^2,3 and ^42,4 respectively. Since the 
Landau- Lifshitz equation is not Galilei invariant, we cannot change the value of v 
by a group transformation. The transformation (3.8) leads to (3.10) with 

K 1? 
= -— + AS'^ + 2{A + ^2) _ . 2^;^ + 2{-B + a) 

X2 = + 45' + 2(^ + 62) + ^ _ + 2{B - a) ^^"^^^ 
Ks = -{2S - vf. 
Eq. (3.6) in this case gives 

«e = y ^ (3-46) 
For S we obtain (3.11) with the constraint 

t/it/2W4- (1-^)' (3.47) 

imposed on these roots. Hence, only even number of roots can be negative (0, 2 
or 4). This however means that all solutions (3.12) . . . (3.20) can occur, though in 
some cases we must impose t/4 < (t/4 = is allowed for v = 2S). 

For S — 0, all solutions (3.22) , . . ., (3.30) can occur. 

4. Algebra ^2,5 

The reduction formula is 

W = R{0 expz[a(0 + ^0 + Bt], ^ = \np + ^(j) (3.48) 

and we must set A = in the Landau-Lifshitz equation. The phase a{^) satisfies 

M + B^ a 



(see (3.6) ). 



18 



The function U = —R^[^) satisfies (3.10) and for 5 7^ we have 

ii'i = i^2 = -f + 452-^^2?^, K^^-AS\ (3.50) 

The solutions that can occur in this case are (3.14) , . . ., (3.20). However, there 
are constraints between various parameters of the solution which follow from the 
requirement of singlevaluedness of W . 

For a = eq. (3.11) again reduces to (3.41) and we only obtain the elementary 
periodic solutions (3.16) . 

For 5" = 0, 7^ we obtain the equation 
Ul^KiU{U -Ui){U 

2a2 (3.51) 

For < f^i < we have t/i < t/2 < and solutions (3.22) , (3.23) , (3.26) 

and (3.27) are obtained. 

For > we have < C/i < C/2 but we obtain no real solutions. 

For Ki < -a^/{h^ + 1) we obtain solutions (3.30) . 

Finally, for 5" = 0, = the solution is C/ — —I^ exp (iV^^) and hence 

i? = i?oexp[±^V^(lnp+i0)] (3.52) 

For S — Ki ^ Q equation (3.10) reduces to 

Ul = K,U{U^ + ^U+ll K, = -^-j0^ (3.53) 



Real solutions are obtained only for Ki < 0. More specifically, solutions (3.27) and 
(3.30) can occur for any Ki < 0. Solution (3.26) for Ki in the range —q^^^ < 
Ki < 0, (3.22) for Ki = — ^pq^jp and (3.23) either for a ^ 0, or Ki ^ — (PTTF 
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5. Algebra A2fi 

The reduction formula is 



W = expi[Q;(0 + Bt + a\nx], ^ = ?7 = arctan - = (3.54) 
and a satisfies 

and = satisfies (3.11) with Ki = K2 ^ + AS'^ + 2a? = -4^2 

For a — Q the equation again reduces to (3.41) . 

A real solution is obtained only for K > SS^, namely (3.16) . It is periodic in 
(f) and hence singlevalued when fi is an integer. 

For 5" = we get a real solution for Ki < 0, namely solution (3.25) which in 
this case reduces to 

U{<P) = - tan2 ^^/\K^\{<P - <Po). (3.56) 



This is a singlevalued function whenever is an integer. 

6. Algebras A2J and ^42,8 

The corresponding reductions lead to equations that we cannot decouple with- 
out introducing higher derivatives, so we will not treat them here. 
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4. Solutions of the Nonlinear Diffusion Equation 



4.1. General Procedure 

Let us consider the system (1.7) and (1.8) , the NDLE for short. We impose 
that the solution be invariant under one of the Lie groups generated by the algebra 
in cq. (2.7) . The functions R and Q will then have the form (3.1) with ^ and (3 
as in Section 3 (different for each subalgebra v42,i, . . . , ^2,8)- 

As for the LL equation, algebra ^2,1 must be treated separately. 

For ^2,2, • • • ,^2,8 we always have (V^^ 7^ 0. Eq. (1.7) and (1.8) reduce to 



(4.2) 

The functions /(C), ^(C), ^(C) and g{^) are defined by the relations 

k _ it- H _ (ve,V/3) 
/ (vo^ ' (vo V 



In order to decouple equations (4.1) and (4.2) , we impose a restriction on the 
functions defined in (4.3) namely 

+ /i = (4.4) 

Eq. (4.1) can then be integrated once to give 

n _i_ /?2\2 

<^^=[S ^ ^2 (4.5) 
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where S is an integration constant. We substitute (4.5) into (4.2) and put 

R{0 = V^im, v = J fim- (4.6) 

The equations are decoupled and the one for U{ri) is aheady in a standard 

r [20] 1 

lorm namely 



with 



(4.8) 



2 B 



/2(V02- 

We now make a further restriction, namely, that M and A^, defined in (4.8) are 
constants. Eq. (4.7) then has a first integral K and we obtain the elliptic equation 
(3.10) with 

Ki^-^ + 4:S'^ + M + N 

X2 = -:| + 452^M-iV (4.9) 

= - 4^2 

In many cases we have M = = and the polynomial on the right hand side 
of (3.10) has a double root at t/3 = t/4 = 1. The solution we obtain for 5 7^ is 
(3.16) with [/4 = 1 ie: 



(C/2 - C/i) sin^ //(?7 - 770) + 1 - C/2 (4.10) 
-5^(1- t^2)(l-[/i), C/i < [/ < C/2 < 



For 5* = 0, > we obtain solution (3.25) i.e 



U ^ - t&n'^ \^{rj-rio) (4.11) 
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4.2. Individual reductions 



1. Algebra ^2,1 
We have 



R^R{t), 



Q ^ ao + ax + by 



(4.12) 



where ao is a constant and Rt satisfies: 



R{1 - R^ 



+ + 62) _ BR 



(4.13) 



Equation (4.13) can easily be integrated (differently depending on whether 
{A -\- B -\- -\- lP'){A — B -\- a? -\- vanishes, or not) and we obtain a 
transcendental equation for R{t). 

2. Algebra ^2,2 

The reduction formula is (3.37) . Wc have m = and (4.4) requires 6 = so 
the solutions are static ones. The variable r] and constants involved satisfy 



Since we have 6 = 0, the solutions are static ones. All solutions (3.12) , . . . , (3.30) 
can occur. 

3. Algebra ^2,3 



and again, the solutions are static ones, since the reduction formula is (3.44). 
All solutions of Section 3.2 can occur. 



77 = In p, M = 
A = B = N = b = Q 



(4.14) 



We have 



rj^x, M = 2(62 + A), N^2B, a = 



(4.15) 
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4. Algebra ^2,4 

The reduction formula is (3.44) and we have 



M =Ar = 0, i; 7^ 



The obtained solutions are (4.10) and (4.11) and they are ^-dependent. 



5. Algebra A2.5 

We have eq. (3.48) with 



77 = ^= lnp+^(/), A^B^N^% 



Since we have M > we obtain the solutions (3.16), (3.18), (3.20), (3.22), 
(3.23), (3.26), (3.27) and (3.30). The parameters of these solutions must 
satisfy, however, certain constraints in order for the solutions to be single- 
valued. For K = 2M we obtain 



2|aH 21aJ_ , 

U = -i^gp^e^"^ (4.18) 



a solution that is not single valued. 

6. Algebra ^2,6 

The reduction formula is (3.54) and we have 



n = J-^ = ^^ A^B^a^M^N^O (4.19) 

so the relevant solutions are (4.10) and (4.11) The solutions are static and 
they are singlevalued for or y^K/8 being integers. 
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7. Algebra A2J 
We put 

W^R{Oexpi[a{0 + ^lnt + b<l>] 

where Ei{x) is the exponential integral function. Moreover, we have b ~ a — 
B ^ M ^ N ^ and the relevant solutions are (4.10) and (4.11) 

8. Algebra ^2,3 
We have 

^ (4.21) 



V 



where ^{x) is the probability integral. We have M — N — A = B = a = 
and so we obtain solutions (4.10) and (4.11) 

We see that time-dependent solutions are obtained for the algebras >l2, 1,^2,4, 
^42, 7 and ^2,3- For ^2,5, A2J and ^2,3 the solutions are trigonometric ones. 

The phases a{^) can be calculated by direct integration, since we have 

1 = -^ 

and U is already known. 

Thus for U given by (4.10) we get 



a{v) = ^ ^^^H^ arctan{ ^ _ || cot /.(^ - ^0) } 
- VW^^^)(^^2-^ arctanj W _ |j cot /i(r? - %) } + ao, 



(4.23) 
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while for U given by (4.11) we have 



a^iJ^S cot{ (77 - 770)} + ao. (4.24) 



5. Conclusions 

The Landau-Lifshitz equation (l.l)has received quite a bit of previous at- 
tention, mainly in the context of continuum Heisenberg ferromagnetic spin sys- 

[21] [22] [23] [24] [26] [26] [27] ™ • ; rr- • ; A i t-. r /-i r.\ 

tems . The amsotropy coefhcients A and B of eq. (1.3)were usu- 

ally set equal to zero. Use was made of the fact that eq (l.l)is integrable, at least 
in the one-dimensional case, or in the two-dimensional, spherically symmetric one. 

[23] [25] 

A one-soliton solution has been obtained a radially symmetric one. In our 
variables R, Q of eq. (1.5), (1.6)this solution corresponds to 

R- ^^ + "1 (Kl) 



where ai and 02 are arbitrary real constants. As noted by Lakshamanan and 

[25] 

Porsezian the soliton spreads in time. 

The solution (5.1) is not among the invariant solutions obtained in this article, 

nor can it be obtained from such a solution by applying transformations from the 

[28] 

symmetry group. As often happens , the method of symmetry reduction that 
does not rely on integrability, provides different solutions for integrable equations, 
than the use of Lax pairs, or Backlund transformations. 

We note that eq. (1.1) with ^4 7^ is not integrable. 

In general, we have reduced the LL equation to the ordinary differential equa- 
tion (3.9). We have integrated eq. (3.9) in terms of elliptic functions whenever p, q 
and m are constants. For the algebras ^2,3, . . . , ^2,6 this was always the case. 
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For algebra ^2,2 we obtained eq. (3.10) only for A = 0, B = b. Let us briefly 
consider the case when the anisotropy coeflicient A does not vanish. We then 



return to the original variable i = p = y/ -\- y'^ = exprj and transform eq. (3.9) 
into 

(5.2) 

For a = 0,b = B this is the equation for the fifth Painleve transcendent Py [20] . 
However, for B ^ b eq (5.2) does not have the Painleve property. According to the 
Painleve conjecture eq. (1.1) is hence, in general, not integrable. 

This has not stopped us from obtaining numerous solutions, both in integrable 
and nonintegrable cases. The algebra ^2,2 (cylindrical symmetry) ior A — 0, B — b 
leads to periodic solutions, as discussed in Section 3.3. The periodicity is in the 

radial variable p. The time dependence is restricted to the phase Q, as is seen in 
eq. (3.37). Moreover the time-dependence is entirely due to the presence of the 
external field B (we have b = B) that generates a rotation between the components 
01 and 02 of the original vector 0. 

Some elementary nonperiodic solutions that we can extract from Section 3 are 



^2 ^ Ui - U2SHU2 - Uijlnp/po)^ 



l + S^U2-Uininp/por 



J.2 ^ ^Rpm - ^1) - U2{U4 - Ui)[2Rl + p^^ + R^p-^^ 



with S, Rq, po, Ui constants and 



:i+i?2)2 

pR'^ 



Q^S I ^ ' dp + acf> + Bt (5.5) 



in both cases. 
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For 5* = we have for instance 



with 



R = [-Ui + J^iln-^r']'^' (5.6) 
V --^1 Po 



R = (5.7) 



Q = a(P + Bt + Qo (5.8) 



For A 0, as mentioned above , solutions are obtained in terms of Pv{p)- 
Their time dependence is again given by the term Bt in the phase Q. 

For algebras ^2,3 and ^42,4 we obtain eq. (3. 10) and a multitude of explicit 
solutions for all values of a, 6, ^4 and B . Note that for 

S^O, B^a, A^-b^ (5.9) 

in particular for the one dimensional (6 = 0), static (a = 0) with no external fields 
{A — B = 0), two of the roots in eq. (3. 11) coincide and the equation reduces to 

= -iS^U - l)2[C/2 + (-^ + 2)U - ^] (5.10) 

Eq. (5.10) only allows elementary solutions like (3.12), ...(3.17), not however 
the elliptic function ones. These occur when the fields A and B are such that (5.9) 
is not satisfied. 

To our knowledge, the NLDE (1.2) has not been investigated from the point 
of view of its integr ability and we have no solutions to compare ours to. 

We have derived many explicit exact solutions of both equations. Looking at 
them we note that most of them have infinite energy. They can describe coherent 
phenomena in various solid state and condensed matter applications. 
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Looking first at the solutions of the LL equation we note that some of our 
solutions have finite energy. In particular, this is the case for (3.43). This solution 
is obtained from the familiar static solution describing n solitons "on top of each 
other" . Its time dependence is given by the factor e*^^ which thus describes a 
rotation of this static solution in the 0i, 02 plane with the angular frequency given 
by the anisotropy B. The other solutions of this class correspond to the static 
elliptic solutions discussed in ref [17] again rotated by e*'^*. 

The solutions corresponding to algebras ^2,1, ^2,3 and ^2,4 have infinite ener- 
gies. As such, they describe various waves in the medium (generalisations of plane 
waves). These can for instance be spin waves; the energy per period is finite. 

An interesting class of solutions are those corresponding to algebras ^2,5 and 
^2,6- Given the choice of parameters, these solutions can be of finite energy; 
however, due to their dependence on the variable = arctan ^, they may become 
singular when x and y vanish. They can be used to describe media with defects. 

Most of the comments made above apply also to the solutions of the NLDE. 
The static solutions in both cases are of course the same. When we consider non- 
static solutions, the most interesting from the physical point of view, are perhaps 
solutions corresponding to algebras ^2,1, ^2,7 and ^2,8- All of them have infi- 
nite energies. The solution corresponding to ^2,1 describes a structure shrinking 
towards the origin (or expanding to infinity - depending on the values of the pa- 
rameters). The other solutions describe field configurations evolving in time. They 
can be used in the description of some physical phenomena in condensed matter 
or solid state physics. 

Among the questions that we plan to return to, we mention the study of "par- 

[32] [33] [34] 

tially invariant" solutions of eq. (l.l)and (1.2), and also "conditionally 

invariant" ones'^^'. A study of solutions involving Painleve transcendents is also 
warranted. 
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